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We theoretically investigate the trion-polariton and the effects of a two-dimensional electron gas 
on its single particle properties. Focussing on the trion and exciton transitions, we set up an 
effective model and calculate the optical absorption of the quantum well containing the 2DEG. 
Including the light-matter coupling, we compute the Rabi splitting and polariton lineshapes as a 
function of 2DEG density. The role of finite temperature is investigated. The spatial extent of the 
trion-polariton is also calculated. We find a substantial charge build-up at short distances as long as 
the Rabi frequency does not exceed the trion binding energy. All our calculations take into account 
the Fermi-edge singularity and the Anderson orthogonality catastrophe. 


I. INTRODUCTION 

Embedding a quantum well (QW) inside a semiconduc¬ 
tor microcavity results in the formation of the so-called 
microcavity polaritons; a coherent superpositon of the 
cavity mode and an exciton, when the cavity is tuned to 
an excitonic transition in the QW. Over the last 15 years, 
these bosonic quasiparticles were shown to have many 
nice properties for studying quantum optics and many- 
body physics in integrated photonic structures. For ex¬ 
ample, their photonic component allows for a direct ob¬ 
servation of the polaritons and a large coherence length, 
while the matter component is responsible for polariton 
interactions. Bose-Einstein condensation of polaritons 
and their superfluid properties have been investigated in 
detail both experimentally and theoretically^. 

The presence of charges inside the quantum well 
that could interact with the polaritons, has been put 
aside until recently. From the theoretical side, the 
mixed electronic-polaritonic system was proposed to 
reach superconductivity at higher temperatures, possibly 
even room temperature, thanks to a strongly attractive 
electron-electron interaction mediated by the polaritons^. 

A pioneering experiment with high density modula¬ 
tion doping was performed by Gabbay et al£, where it is 
shown that strong light-matter coupling is possible, even 
though neither an excitonic nor a trionic state could be 
resolved. Recently, in the experiments by Smolka et al.— 
polariton formation was observed over a range of electron 
densities, both with and without a magnetic field. 

For what concerns the effect of the electron gas 
on the polaritons, various contributions both from the 
experimental^ - — and theoretical^ side have been made. 
Those works considered the electron-polariton scattering 
physics, but it was assumed that the polariton quasipar¬ 
ticle itself was not affected by the electron gas. 

The optical excitation of an electron gas (metal or 
doped semiconductor) is an important problem in many- 
body physics that was introduced by Mahan^ and to 
which Nozieres, De Dominicis^ and Andersonii made 
seminal contributions. Nowadays, the Fermi edge singu¬ 
larity is still an active topic of research that has recently 


attracted the attention of the community working on ul¬ 
tracold atomic gasesi2r— . 

The role of the electrons on the properties of the 
exciton-polariton itself for a hi ghly do ped quantum well 
has been investigated in Refs. [15] and [16]. In our previ¬ 
ous worbii, we have discussed the single particle exciton- 
polariton properties, properly taking into account the 
many-body physics involving the Anderson orthogonal¬ 
ity catastrophe (AOC) and Fermi edge singularity (FES). 
Unfortunately, an exact evaluation of the optical response 
was only possible for negligible electron-electron interac¬ 
tions. In the low density regime, where trion correlations 
dominate, this makes it a poor approximation. In the 
case of a low density electron gas, trion-polaritons have 
been observed and discussed theoreticall y 4 : 18 ! 19 . A theo¬ 
retical model for trion-polaritons, that takes into account 
the AOC and FES, is however still missing. It is the pur¬ 
pose of this paper to fill this gap. We will focus on the 
linear properties of trion-polaritons, deferring the study 
of polariton-polariton interactions to a future work. 

This paper is organized as follows: in Section [TT] the 
notations used in this article are introduced and we set up 
our model and show how to calculate the optical response 
of the two-dimensional electron gas. Section Hill presents 
some linear polariton properties such as the polariton 
energies, lineshapes. In the same Section we calculate 
the spatial distribution of electrons in the lower polariton 
state. The last section HVl gives a summary and outlook. 


II. HAMILTONIAN AND OPTICAL RESPONSE 
OF THE QUANTUM WELL 

A. Model set-up and formalism 

The system under consideration is a quantum well con¬ 
taining a non-interacting two dimensional electron gas 
(2DEG), placed inside a planar microcavity. For the sake 
of simplicity we consider a spin polarized 2DEG (the gen¬ 
eralization of our results for an unpolarized electron gas 
is straightforward). When a photon impinges onto an 
empty QW, an interband electron-hole pair can be ere- 
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ated and forms a bound electron-hole pair, the exciton. 
In the presence of a 2DEG, the exciton will interact with 
the electrons, with the possibility of forming a trionic 
bound state, where the valence band hole is bound to 
two electrons. Typically, the trion binding energy is much 
weaker than the exciton binding energy. Therefore, the 
trion can be approximated as an electron that is bound 
to an exciton through an effective clectron-exciton poten¬ 
tial V x ~ e . This approximation reduces the trion physics 
from a three-body to a two-body problem. Within this 
approximation, it becomes feasible to treat the interac¬ 
tion of the exciton with all the electrons in the 2DEG. 
Provided that the electron-electron interactions can be 
neglected and the exciton mass tends to infinity, the for¬ 
malism of Combescot and Nozieres^S yields an exact de¬ 
scription of the optical response properties. We will re¬ 
strict in this paper to these approximations and leave 
the inclusion of electron-electron interactions for future 
work. 

We thus start from the following Hamiltonian: 


plane waves \k)). Using the Heisenberg picture for the 
time dependence of the operators we can write it as 

G(t) = (FS|e ijfft ^ A -e“ iff ^t|FS) 

= (FSIe-^-^lFS). (3) 


The second line expresses that the time evolution of the 
electrons after the injection of the exciton is governed 
by the modified Hamiltonian H = Hm($x'$x = 1), for 
which |FS) is no longer an eigenstate. The overlap with 
the unperturbed Fermi sea will therefore become time- 
dependent, and in particular for long times it will decay 
as a powerlaw. The value Eq = ^2 k<kp e k is the ground 
state energy of |FS). The expectation value ([3]) can be 
straightforwardly computed numerically 12 ' 20,21 - as 


G(t) = det 


1 — flF + flF^(t) 


( 4 ) 


Here, uf is the Fermi-Dirac distribution. The matrix 
A (t) is given by 


H = H m + Vlm , 


with 


hk'(t) = Y,m{p\k')e- i(e »- ek)t , ( 5 ) 

p 


H m = exi’x'&x + 55 + 55 V kk’ e Ck'C^x$x- 

k k,k' 


where one has H\p) = e p \p). Finally the 2DEG absorp¬ 
tion is given by 


Here, Hm describes the matter degrees of freedom of the 
quantum well. Because the exciton (created by '(’x) is 
assumed to have infinite mass, it is sufficient to consider 
a single exciton mode that is localized in real space with 
energy ex- The kinetic energy for the electrons, created 
with 4, is given by e k = k 2 (units are h = 1, m e = 1/2). 
The last term in Hm describes the electron scattering 
by the exciton. In particular we will use an attrac¬ 
tive exciton-electron potential V x ~ e < 0, for which it 
is known that in two dimensions this always results in 
the presence of a bound state with energy £t < 0. Fur¬ 
thermore, we restrict ourselves to the spheric symmetric 
l = 0 angular momentum channel. 

The coupling of the QW to the photon field is treated 
semiclassically as 

V LM = gA L e~ iULt -ipx + h.c. (1) 

Here g is the coupling constant between the optical mode 
in the cavity and the exciton center-off-mass at zero in¬ 
plane momentum k = 0. The amplitude Al represents a 
coherent drive by an external laser field. 

In order to calculate the optical response function of 
the QW, we start from the free 2DEG and treat the light- 
matter coupling term as a perturbation to the system. 
The optical response G{t) of the QW is then given by 
linear response theory as2£ 

G{t) = (FS|T{<Mt)4}|FS> (2) 


i r°° 

A{u) = - Re d te-* ut G(t). (6) 

n Jo 

It has been shown^S that in the limit t ej, 1 the 
function G(t) is given as a sum of two powerlaws, 

giu)2t 

g(t»e-/)=C 1 — + c 2 — (7) 

with Ci t 2 some constants. The powerlaw decay is a 
manifestation of the so-called Anderson orthogonality 
catastrophe^ which states that the ground state of the 
system with and without the scattering potential (due to 
the exciton) are orthogonal to each other. The powers 
0 : 1,2 can be related to the scattering phase shift at the 
Fermi level of the conduction electrons scattering off the 
potential: 



(SfV 

Oi = 

— , 


\ir J 


(S F / 

02 = 



( 8 ) 


The phase shifts satisfy Sf £ [0,7r], making 01,2 € [0,1]. 
Furthermore, since the Laplace transform of the separate 
terms in 0 is given by 


/•OO 


/ d te~ iujt 

c — 

Jo 

t a _ 


|CjT(l - a)e- i! r 
{u — n + ip + ) 1 ~ a ’ 


( 9 ) 


where T is the time-ordering operator, and |FS) is the with T(a;) the gamma function, it follows that the absorp- 
unperturbed Fermi sea (a Slater determinant built with tion spectrum contains two singularities, the so-called 
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Fermi edge singularities. The positions of the singulari¬ 
ties in frequency domain (‘thresholds’) is given by2£ 

0Ji=e T + A (e F ), 

0J2 = A(£f) + e_F, (10) 

where A(e^) = Y2 P< k F ( £ p — £ p ) is the energy difference 
between the system’s ground state before and after the 
quench. Here, the single particle energies e p are obtained 
by diagonalizing H. 

Let us now be more specific about the choice for 
the exciton-electron potential: in the remaining part of 
this paper we work with an attractive Yukawa potential 
V x ~ e (r) = ae~ r l b /r with a < 0,b > 0 numerical con¬ 
stants. The thresholds (fTUl) and the corresponding pow- 
erlaw exponents in real time domain (©) for the above 
choice of V are depicted in Fig. |T] as a function of Fermi 
energy. Because the scattering potential is attractive, the 
spectrum e n contains a bound state with energy et < 0, 
corresponding to the state in which an electron is bound 
to the exciton, i.e. a trion. The trion binding energy 
is thus given by the amount of energy needed to disso¬ 
ciate the trion into an exciton and an electron. In the 
following, the binding energy will be used as an energy 
scale. 

At zero density, the lowest threshold corresponds to 
the trion and the upper threshold is given by the exci¬ 
ton energy. For sake of simplicity, in the following we 
will call the lowest threshold the trion, and the high¬ 
est threshold the exciton for all densities. Still at zero 
density, we find that the powerlaw exponent in time do¬ 
main corresponding to the trion (op) is exactly one, while 
the exciton exponent 0:2 is zero. From eq. © we then 
immediately see that the exciton corresponds to a delta 
function in absorption while the trion is completely flat 
in frequency domain. For increasing electron density, the 
trion exponent gets larger than the excitonic exponent, 
the crossover being for the Fermi energy comparable to 
the trion binding energy. The latter means that the trion 
tends to get more delta-function like in absorption for in¬ 
creasing 2DEG density, while the exciton gets flattened 
out. This behaviour has been numerically verified to hold 
for several short-range potentials. In particular, we also 
checked it for the more realistic 1/r 4 exciton-electron po¬ 
tential. The two thresholds are getting more separated 
for increasing Fermi energy. This would correspond to 
an increase of the trion binding energy for higher elec¬ 
tron densities. This seems unphysical and should be at¬ 
tributed to the neglect of electron-electron interactions. 

Note that, in order for our model (the exciton as ele¬ 
mentary boson) to be valid we must avoid very high elec¬ 
tron densities. Specifically, we are restricted to densities 
satisfying ep -C \ex\- This means that the electron inter¬ 
particle distance should always be larger than the exciton 
Bohr radius. Only in this situation it is not possible for 
the electrons to see the internal structure of the exciton. 
Typically, in GaAs QW’s it holds that et — O.lex, so in 
order for our model to stay valid we must have approx- 




Figure 1. a) Thresholds wi ,2 as a function of Fermi energy, b) 
Powerlaw exponents 01,2 as a function of Fermi energy. Both 
exponents are seen to be smaller than one, corresponding to 
two singularities in the absorption spectrum. For Fermi en¬ 
ergies larger or comparable to the trion binding energy, the 
trion exponent (blue open circles) starts getting smaller than 
the exciton exponent (red squares). 


imately Ef < 10|£t|- This explains the horizontal range 
on the figure© 


B. Lineshapes, oscillator strength, temperature 

For the Yukawa potential mentioned earlier, we have 
numerically computed the response function G[t) in 
the time domain up to times for which the asymptotic 
regime, eq. 0 - sets in. In Figs, ©a-c) the modulus 
of the response function is plotted for increasing 2DEG 
density from top to bottom. Note the double logarit- 
mic scale to evidence the powerlaw nature of the decay. 
Solid lines correspond to the numerical simulation while 
the dashed (dotted) line correponds to the powerlaw de¬ 
cay using the trion (exciton) exponent from Fig. [lib). 
To obtain the absorption, we numerically performed the 
Laplace transform of the time series using an appropri¬ 
ate damping term to avoid manifestations of the Gibbs 
phenomenon. The corresponding absorption spectra are 
seen in Figs. ©d-f). All spectra are shifted with respect 
to the trion threshold u>i. The spectral lines are seen 
to have an asymmetric lineshape with a powerlaw tail 
at the high frequency side of the threshold. This is of 
course due to the powerlaw exponent in eq. ©. Only 
for the real time exponent being exactly zero, the line- 
shape becomes a symmetric Lorentzian. In theory, the 













4 




10° 


=3 

Hi,. _Q 



8 p /8 T =1 


V/ VI 





0 

v |! 

im 1 




.A' 2 

(b) 

Q. 

O 

GO 



^ (e) 


10° 10 1 10 2 < 0 2 4 



8 f / 8 t =3 


(f) 

JL_ 

0 2 4 

(co-co 1 )/|e T | 


Figure 2. a-c) Solid lines: modulus of the response function 
|G(t)| on a double logaritmic scale for increasing Fermi en¬ 
ergy. The dashed (dotted) lines depict the powerlaw decay 
using the exponents a\ ( 02 ) from eq. lf8)l. (d-f) Absorp¬ 
tion spectra obtained after performing the Laplace transform 
of the time series. Spectral weight is shifted from the exci- 
ton singularity at low densities towards the trion threshold at 
higher densities. 


absorption on the left hand side of the trion singularity 
should be identically zero. Due to the numerical damp¬ 
ing, this sudden step is slightly rounded. Furthermore, 
for increasing Fermi energy we see that the distance be¬ 
tween the two thresholds indeed increases, as mentioned 
before. The trion is seen to get narrower again for higher 
densities (the corresponding real time exponent goes to 
zero for high densities, see Fig. nh). Experimentally 
many effects such as electron-electron interactions, finite 
temperature, will nevertheless broaden the spectral line. 

To quantify which of both singularities is most dom¬ 
inant in the spectrum, we define the oscillator strength 
as the integral of the absorption spectrum around each 
singularity. The boundary between the two thresholds 
is taken as the frequency where the absorption has its 
minimal value in between. Since from the definition © 
we have JM° dujA(uj) = 1, each oscillator strength is less 
than 1. Fig. [3] shows the oscillator strength as a func¬ 
tion of Fermi energy. At zero density all spectral weight 
is in the exciton (red squares), that is a delta function. 
This is indeed the only transition that is possible in our 
model, since at least one electron is needed to form a 
trion. For increasing Fermi energy, oscillator strength 
is transferred from the exciton towards the trion (blue 



Figure 3. Oscillator strength of the trion and exciton as a 
function of Fermi energy. At zero density the exciton contains 
all oscillator strength. The trion starts dominating for Fermi 
energies larger than the trion binding energy. 


circles). Both singularities have equal oscillator strength 
when the Fermi energy is comparable to the trion binding 
energy. For very high 2DEG densities most oscillator is 
gathered in the trion. This is consistent with the spectral 
lines in Figs, ©d-f) where the exciton is seen to be dom¬ 
inant at low densities, while the trion starts dominating 
for high electron densities^. 


Because of the simple form of eq. ©, it is straight¬ 
forward to investigate finite temperature effects on the 
spectral lineshapes. In the Fig. ©1 we show the time evo¬ 
lution of G(t) for a fixed Fermi energy ef/et = 1, and 
for several temperatures T/et■ For finite temperature, 
the time series are always seen to coincide with the zero 
temperature result for short times. However, deviations 
from the T = 0 start showing up at the time proportional 
to the inverse temperature, i.e. t ~ T _1 . For long times 
the finite temperature series decay exponentially to zero 
(not visible on the double logaritmic plot). Figure ©3 
depicts the corresponding absorption spectra. For T = 0 
we have the same lineshape as in Fig. ©. For increasing 
temperature, the spectral line shapes become broadened 
and tend to get more Lorentzian. As expected, increas¬ 
ing the temperature too much, causes the two thresholds 
to merge. 
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Figure 4. a) Time dependence of |G(t)| for a fixed density 
ef/et = 1 and for several temperatures T/et = 0, 0.1, 1, 5. 
For finite T, the response function exponentially decays, with 
a decay time set by the inverse temperature, b) Absorption 
spectra obtained by Laplace transformation of the time se¬ 
ries. Finite temperature causes a broadening of the absorp¬ 
tion peaks and makes them more symmetric, i..e Lorentzian. 


III. TRION-POLARITON 


A. Photon spectral function 


We can now turn on the light-matter interaction be¬ 
tween the cavity photon and the charged quantum well. 
In particular, we are intrested in the linear properties of 
the coupled light-matter system, i.e. we do not consider 
photon-photon interactions. The main quantity of intrest 
is the photon spectral function, given by 


D{u) = — Im- 

w — 


1 

W c + i\g LM \ 2 G(iuj)' 


Here uj c is the cavity mode energy and G{iu >) the Laplace 
transform of eq. ©• The microcavity is assumed to be 
lossless and we only consider a single photon mode at 
normal incidence. 

From now on, we restrict our attention to the trion 
in the 2DEG absorption spectra at zero temperature. In 
particular, as the trion energy shifts with increasing elec¬ 
tron density (see Fig. [lji), we can put the photon into 
resonance with the trion for every density. The result¬ 
ing lower polariton energy Epp is seen in Figs. [5ji, [6ji 
for two different values of the light-matter coupling, i.e. 

= 0.5|£t| vs <7lm = 3|£t|- In both cases, we see a 


monotonic increase in the trion-polariton ‘Rabi splitting’ 
as a function of Fermi energy. This can be attributed 
to the gain of trion oscillator strength as the Fermi en¬ 
ergy increases. For low electron density, the increase of 
the trion oscillator strength can be simply understood 
in terms of an increased trion-photon overlap^. In re¬ 
cent experiments^, a decrease of the Rabi splitting was 
however observed for higher electron density. Introduc¬ 
ing hole recoil and electron-electron interactions could 
solve this discrepancy. In the subfigures (1-3) the po¬ 
lariton lineshapes (red solid) are depicted for some fixed 
Fermi energies, the values being indicated by the black 
dashed lines in the large figures. For comparison, the 
2DEG absorption has also been depicted in blue dashed 
lines. For small light-matter coupling, the exciton does 
not influence the cavity mode. The lower polariton now 
truly corresponds to the trion-polariton as being a su¬ 
perposition of trion and photon. Note that in principle 
the lower polariton should have zero linewidth, but for 
visibility we gave it some small value. The other polari¬ 
ton modes intrinsically possess a finite linewidth due to 
the finite absorption of the 2DEG at the corresponding 
polariton energies. For larger light-matter coupling, the 
photon starts admixing with the exciton resulting in a 
non-zero quasiparticle shift at the lowest electron density 
under consideration. We also see a clearly visible polari¬ 
ton mode above the exciton threshold. Furthermore, in 
Fig. [61(2) the middle polariton is seen to have an asym¬ 
metric lineshape with a large tail towards the low fre¬ 
quencies. This feature stems from the interplay between 
the two skewed lineshapes from both the trion and exci¬ 
ton and the fact that they both contribute equally to the 
polariton formation because of their oscillator strengths 
being the same at the corresponding Fermi energy. 


B. Electron density in the trion-polariton state 


One of the attractive features of microcavity polari- 
tons is their strong optical nonlinearity, allowing the cre¬ 
ation of an interacting polariton quantum ‘fluid’. A Fes- 
hbach resonance mechanism for enhancing the polariton 
interactions^ has recently been investigated^. Unfortu¬ 
nately, in the current state of the art microcavities, the 
single photon nonlinearities are too weak in order to enter 
the photon blockade regime. 

It has been recently suggested that the electrons in the 
QW could enhance interactions between the polaritonsi 
In the previous part of the paper, we have shown the 
lower trion-polariton to be the good quasi particle for a 
highly doped QW embedded in a planar microcavity. Be¬ 
cause the matter component of the polaritons is respon¬ 
sible for the interactions, it is instructive to investigate 
the spatial structure of the trion-polariton. In particular, 
we will calculate the electron density in the lower polari¬ 
ton state through the electron-exciton density correlation 
function: 
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Figure 5. a) Lower polariton energy as a function of Fermi en¬ 
ergy, obtained by putting the cavity mode into resonance with 
the trion for all Fermi energies. The light-matter coupling was 
taken qlm = 0.5|£t|- (1-3) Polariton lineshapes (red solid) for 
some fixed Fermi energies (black dashed lines in Fig. a). For 
reference, the 2DEG absorption is also shown in blue dashed 
lines. The exciton threshold is not visible because it lies at 
higher energies. 



Figure 6. Same quantities as in Fig. [5] but now for a light- 
matter coupling qlm = 3|et|. The larger light-matter coupling 
now admixes the photon with both the trion and the exciton, 
resulting in an upper polariton with an energy higher than the 
exciton energy. In (2) the middle polariton is seen to have an 
asymmetric lineshape, inherited from the 2DEG absorption. 


_ (h e (r)nx (0)) 

<M0)> ’ 


(ii) 


with n(x) = (x)ip(x) the density operator at position 
x. Because the exciton is assumed to be infinitely heavy, 
we can place it in the origin without loss of generality. 
Since we only consider s-wave scattering, the meaning of 


the spatial coordinate r is the relative radial distance 
from the electrons to the exciton center-of-mass. We 
start again from the unperturbed Fermi sea and consider 
the Hamiltonian Vlm given by Eq. 0 as a small pertur¬ 
bation to the 2DEG. Using Kubo’s response theory, we 
find that the lowest order contribution to the above cor¬ 
relation function is of second order in the perturbation, 
i.e. proportional to |gHi| 2 . It is given by 


(n e (r)nx(0)) = \gAA 2 lim lim [ dse r,s [ dr e r,T e lulL ^ s (FS|'i/>x(0, s) [n e (r, t)nx(0, t)l V’tUO, r)|FS) + h.c. 
ri^0+to^-ooj t J t 


/ oo 1 

dug nm {u) - - — 

-oo W - (wi 


- (w L -e n )oJ- ( UJ L - £ m ) ' 


( 12 ) 


Here, we adiabatically turned on the laser in order to 
find the stationary state of the electron density under the 
continous wave excitation (see details in the appendix). 
The wave functions </> n (r) are the single particle eigen¬ 
states of the Hamiltonian H in position space, satisfying 


Hcj) n (r ) = £ n <p n (r). The matrix g ra m(w) is given as 
gnmijx') — 

7r 

9nm{t) - ^C(fc)<)>m(<?)(FS|c+e- l ( fi - Bo ) t c fe |FS)(13) 

k,q 


Re 


Ate lujt g n 


*(*), 
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The expectation value in the above equation looks similar 
to eq. (0 and it can again be calculated numerically^. 
This means that our results on the electron density prop¬ 
erly take into account the Anderson orthogonality catas¬ 
trophe and the Fermi edge singularity physics. 

The expectation value (fill) should be calculated in the 
lower polariton state. In the resulting expression eq. m 
the polariton shows up through the laser frequency u^. 
In particular, we consider resonant continous wave ex¬ 
citation, thus with the frequency of the laser resonant 
with the lower polariton energy £lp- Note that erp is 
a free parameter in this calculation, but its actual value 
has been calculated in the previous part of this paper, 
see Figs. 0i, 0i. 

The exciton density in the lower polariton state can be 
calculated analogously. It is found to be 

<M0)> = M 2 r du A{u}) (14) 

J -oo (w - Wl ) 2 

where the absorption A(u>) is given in eq. 0. 

The electron density along the radial direction and in 
the lower polariton state is shown in Fig. [7] for several 
Fermi energies and polariton energies. The trion Bohr 
radius ap is taken as the lengthscale corresponding to 
the bound state ~ exp(—r/ap) with energy £p. 

For the lowest Fermi energy, a small Rabi frequency 
(blue solid lines), defined as the energy difference between 
the lower polariton and the trion threshold uq, is seen 
to increase the electron density at short distances from 
the exciton. This can be understood as follows: when 
the photon is annihilated and an exciton is created, the 
Fermi sea is shaken up due to the appearance of the ex¬ 
citon. For small Rabi frequency, it takes longer time for 
the exciton to be reconverted into a photon. This means 
that there is a longer time for the electrons in the Fermi 
sea to adjust themselves to the presence of the exciton: 
a screening cloud of electrons is able to build up near the 
exciton. The reason for the increase of electron density 
is the attractive exciton-electron potential, i.e. the pres¬ 
ence of the bound state. For larger Rabi frequency, there 
will be much less time for the electrons to reorganize; as 
soon as they start readjusting, the exciton has already 
dissapeared and turned into a photon. This is seen by 
the green dashed-dotted line where the exciton density 
tends to flatten. For large Fermi energies (comparable 
to or larger than the trion binding energy), the electron 
density does not depend too much on the Rabi frequency. 
We argue that the time it takes for the individual elec¬ 
trons to adjust to the presence of the scattering potential, 
is proportional to the inverse Fermi energy. The larger 
the Fermi energy, the less time it takes to reorganize the 
electrons at the Fermi level and the system will always 
be able to build up a screening cloud. Nevertheless, the 
absolute value of n e (r)/no decreases for increasing Fermi 
energy, simply because there is less time to screen the ex¬ 
citon. In all cases, the electron density oscillates and the 
envelope decays as 1/r 2 at large distances, corresponding 
to the Friedcl oscillations. 
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Figure 7. Electron density in the lower polariton state (nor¬ 
malized with respect to the unperturbed density no = kp /47r) 
as a function of radial distance to the exciton center-of-mass. 
At large distances, the density stays unperturbed as the cor¬ 
relation function always asymptotically goes to 1. For a high 
density 2DEG the electron density does not depend too much 
on the Rabi frequency because the electrons screen the exci¬ 
ton on a time scale of the inverse Fermi energy. For smaller 
Fermi energy, a low Rabi frequency still allows the electrons 
to screen the exciton. For higher Rabi frequency, the elec¬ 
trons do not have enough time to react onto the presence of 
the exciton so that the electron density is unperturbed. 


From the above, we see that there should be a fre¬ 
quency window in which the trion gets screened on a 
length scale comparable to the bare trion Bohr radius, 
making it a truly localized, electrically charged particle, 
while still ensuring strong coupling. The Rabi frequency 
should be small enough such that there is time for the 
electrons to form for the trion, by binding to the exci¬ 
ton. Since the formation time typically goes as 1 /et, 
one should maintain Qr < et- On the other hand, one 
should avoid running into the weak coupling regime by 
taking the Rabi frequency too small. 


IV. CONCLUSION AND OUTLOOK 

We have set up an effective model describing the for¬ 
mation of trion-polaritons in a microcavity embedding 
a doped quantum well. Taking into account impor¬ 
tant many-body effects such as the Anderson orthogo¬ 
nality catastrophe and Fermi edge singularity, we were 













able to qualitatively describe the 2DEG absorption, the 
crossover of oscillator strength from exciton to trion and 
the polariton lineshapes. The role of finite temperature 
has been discussed as well. Finally, we calculated the 
electron density profile in the lower polariton state and 
concluded that for a specific range of Rabi frequencies, 
the trion-polariton behaves as an electrically charged lo¬ 
calized particle. 

A major simplification in our model was the neglect 
of electron-electron interactions. Including this physics 
would probably solve the decreasing linewidth of the 
trion for increasing electron density. Also, the binding en¬ 
ergy of the trion then becomes a function of Fermi energy 


instead of being a fixed parameter. Most importantly, a 
self-consistent treatment of electron-electron interactions 
in the calculation of the electron density profile is needed 
to see whether a charge build-up near the exciton is still 
possible. 
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Appendix A: Exciton-electron density correlation function 

In this appendix we derive, using Kubo’s reponse theory, the expression for the electron-exciton density correlation 
function. 

Starting from the unperturbed Fermi sea and treating the Hamiltonian dTj) as a perturbation to the system, we need to 
go second order in this perturbation to have a first non-zero contribution to the density-density correlation function. 
Using standard Kubo response theory we find 


{fie (r,t)h x (0,f)) = \gA L \ 2 


ft 

f s r 

/ ds 

/ dr 

' to 

'to L 


o^l(s T )(FS| %j)x (0,s) he ( r , t) fix(0,t) V’x (0,r) |FS) + h.c. (Al) 


where the remaining expectation values need to be calculated in the unperturbed state. 

We perform the calculation in two dimensions. Restriction to the spheric symmetric states comes down to using the 
proper Hankel transform pairs. Using the interaction picture for the time evolution of the operators we obtain 

(FS| i>x (0,s) fie (r, t) nx{ 0, t) & (0,r) |FS) = e iE ^- T \FS\e~ iRs ^ e (r, t ) (r, t) e li?T |FS), 

where H is the Hamiltonian governing the dynamics of the electrons in presence of the exciton scattering potential 
and Eq is the ground state energy of the Fermi sea without the exciton. 
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Using the Fourier transform for the system enclosed in a box of area S , we obtain 

(FS|e- l5 >t ( r; t ) ^ ( r , t ) e ^|F S ) = ^Y (FS| e - lfis 4 (t) c q (i) e l5T |FS), (A2) 

* k,q 

where the plane wave creation operators satisfy {ck, c^} = <5k.q- The time evolution of the operators is governed by 
the Hamiltonian in presence of the exciton: 

H = ^^£fccj.Cfc + Y^ V kk , c\,c k . (A3) 

k k,k' 

The eigenstates are no longer plane waves and we therefore make a unitary transformation from the plane wave basis 
to the basis of eigenstates of H , we call them scattering states. The latter are characterized by their principal quantum 
number n and the angular momentum l , because the scattering potential under consideration is spheric symmetric. 
This gives 


i Y e i(k - q)r (FS|e-^ s 4 (*) c q W e i/?T |FS) (A4) 

* k.q 

1 oo oo 

= ^E ei(k_q)r E E (ci\n,l)(m,p\k}(FS\e-^cl p (t)c n i(t)e^\FS). (A5) 

k.q n,m— 0 l,p——oo 

Since the Hamiltonian H is diagonal in the operators cL we have a trivial time evolution. We than have 


OO OO 


_£ e i(k- q)r £ Y (^\n,l)(m,p\k)(FS\e-^cl p (t)c nl (t)e inT \FS) 


k,q 


n,m=0 l,p——oo 


= ^E e ‘ (k " q)r E e iSmt e~ iSnt Y (q|n,/)(m,p|k)(FS|e-^ct ip c ni e^|FS) 


(A6) 

(A7) 


k,q 


n,m =0 


l,p=— oo 


= ^E ei(k ” q)r E e igmt e~ iSnt Y {q\nJ)(m,p\k)(FS\e- iH f s cl p e iSs e- iS ^- T '>e- iST c nl e iST \FS) (A8) 


k,q 


n,m=0 


l,p——oo 


= ^E el(k " q)r E (? imt <r iSnt Y (<l\nJ)(m,p\k)(FS\cl p (-s)e- iB ^c nl (-T)\FS) 


k,q 


n,m =0 


l,p=— oo 


= sE el(k ' q)r E e igmt e~ ignt e~ iSmS e iSnT Y faK l) <m,p|k) (FS^e"^'-^^), 


(A9) 

(A10) 


k,q 


n,m =0 


l,p—— oo 


where in the third line we introduced the unit operator. The initial state |FS) is built with plane waves but the 
operator c n i destroys a scattering state. We therefore make again a transformation from scattering states to plane 
wave states, finally yielding 


(FS\e~ iSa Tpi (r ,t) iPe (r,t)e li?T |FS) 

1 OO 

= - E e i(k " q)r Y e igmt e~ iSnt e~ igm 

k.q n,m =0 


s e ie n r 


(All) 


E {^\n,l)(m,p\k) Y (n,/|K)(Q|m,p)(FS|c^e-^ s - T )c K |FS). 

l,p=— oo K.Q 


The sum over k.q defines the wave function in position space. Writing (k|n,/) = cj) n i( k), the eigenstates of the 
Hamiltonian H in plane wave basis, we define 


</>ni(r) = ^E e * kr ^ni(k)) 
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and we have 

(FS|V>x (0,s)n e (r,t)hx{0,t) ^ (0,r) |FS) (A12) 

= e iEoi8 - T) (FS\e~ iBs i>l (r, t) iji e (r, t) e lflr |FS) (A13) 

oo oo 

= E E ^ n| (r)^mp(r) e is mt e -is n t e -ie m s e iE n T ^* ; (K)</> m p(Q)(FS|cqe _i ^ _ ' Eo ^ s_r ^c K |FS) (A14) 

n,m—0l,p——oo K.Q 

oo oo 

= E E (A15) 

n,m —0 l,p——oo 

Here, we introduced shorthand notation for the double sum over K, Q (dropping the angular momentum indices in 
order not to overload the notation): 


9nm(t) = E ^(K)</» mp (Q)(FS|c^e- i ^- B “) t c K |FS). 

K,Q 


(A16) 


This expression is exactly the one from eq. (1131) in the paper. The above expectation value should be plugged into 
equation dm This gives 

(n e (r,i) n x (0,i)) (A17) 

OO OO pt pS 

= \gA L \ 2 E E ^«C(r)e^e^‘/ da dr g nm (s - t) + h.c. (A18) 

n,m—0 l,p= — oo to •* to 

The hermitian conjugate term can be shown to be the same as interchanging s and r, yielding 
(h e (r ,t) h x (0 ,t)) 


(A19) 


OO OO pi pZ 

= \gA L f E E ^(r)C P (r)e fe " mt e- fe '“‘ ds dr e^ s ~A g nm (s - r), (A20) 

n,m—0l 1 p— — ac ^ to 

where now we can both have s < r and s > r. Introducing the Fourier transformation of g nm , we can now easily do 
the integration over s, r. Furthermore, in order to obtain the stationary state under the continous excitation, we add 
a small exponential term, meaning we slowly (for example, as compared to the lifetime of the cavity mode) turn on 
the laser. If we then take to —> — oo, this will eliminate all transient behaviour, originating from the initial time to- 
In other words, we do not want our result to depend on to- 
So, introducing the Fourier transform for a function p(f), 


p(t)= / due~ lujt g{co) 

J — OO 

i r°° 

B M = ^ J dt e lult p (t), 


we obtain, by adiabatically turning on the laser (only focus on the time integrals) 


pt pt pOO 

e ig m t e -i£ n t Hm Hm / dse r,+ s / dTe y+r e iuU S -T) e -ie m s e ie n T / e ~M‘-r) ( w ) 

^0+J tQ J tn J_ QO 


t pt 

■ + c 


= e iSmt e- ignt lim lim / dujg nrn (u})l dse^ s [ dr e n+T e ^ (s - T ) e - iE ~ m V £ ~" T e ~ iul{s - T) 

V—>0+ to - >■—oo J _ QO J tQ J to 


diO g n m (u) ' 


1 


- (w L ~ Em) W - (wl - £ n ) ' 


The last line follows from straightforward integration and taking the proper limits. Restoring the double sum over 
n, m, we obtain equation m in the article. The exciton density can be calculated completely analogously. 



